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Abstract 



We give a new method for the evaluation of a class of integrals of rational symmetric functions 
in N pairs of variables {{xa,ya)}a=i,...N arising in coupled matrix models, valid for a broad class 
of two- variable measures. The result is expressed as the determinant of a matrix whose entries 
consist of the associated biorthogonal polynomials, their Hilbert transforms, evaluated at the 
zeros and poles of the integrand, and bilinear expressions in these. The method is elementary 
and direct, using only standard determinantal identities, partial fraction expansions and the 
property of biorthogonality. The corresponding result for one-matrix models and integrals of 
rational symmetric functions in N variables {xa}a=i,...N is also rederived in a simple way using 
this method. 
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1 Integrals of symmetric rational functions 



Let dfj,{x,y) be a measure (in general, complex), supported on a finite set of products of 
curves in the complex x and y planes, for which the semi-infinite matrix of bi-moments 
is finite: 

:= y"d/i(x,y)a;V < oo, 0, Vj, A; G N. (1.1) 

The integrals are understood to be evaluated on a specified linear combination of products 
of the support curves. Assuming that, for all > 1, the N x N submatrix {Bjk)o<j,k,<N-i 
is nonsingular, the Gram-Schmidt process may be used to construct an infinite sequence 
of pairs of biorthogonal polynomials {Pj{x), Sj{y)}j=o...oo, unique up to signs, satisfying 



diJ,{x,y)Pj{x)Sk{y) = 6jk, (1.2) 
and normalized to have leading coefficients that are equal: 

P,(x) = ^ + 0{x^-'), S,{x) = ^ + Oiy^-'). (1.3) 
We will also assume that the Hilbert transforms of these biorthogonal polynomials, 

PnifJ') ■= [ dij,{x,y) ^"'^^\ Sniv) ■■= [ dij,{x,y) ^"'^^\ (1.4) 

exist for all n G N. 
For > 1, let 

r r ^ 

:= / dfi{xuyi) ... dfi{xN, yN)AM{x)A{y) = A^! JJ /i„ (1.5) 

^ ^ n=0 



where 



N N 



An{x) ■=Ylixi- Xj), AN{y) ■=Yliyi- Vj) (1-6) 

i>j i>j 

are Vandermonde determinants. Such integrals are of particular interest in two-matrix 
models IH 13 El UHl HI] and may be interpreted as the reduction to the space of eigen- 
values of the integral defining the partition function on an ensemble of pairs of N x N 
matrices having a probability measure that is invariant with respect to conjugation by 
unitary matrices, and admitting a Harish-Chandra-Itzykson-Zuber type reduction 
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to the 2N dimensional space of eigenvalues. The joint probability distribution (in the 
case of a real Borel measure) is given, after normalization, by the integrand of ()1.5j) with 
eigenvalues (xi, . . .xn), {yi, . . .i/n) having values on the support curves of the measure 
dii{x,y). 

Now choose four sets of complex constants {^a, Cdi Vi^ /ifc}a=i,...ii,/3=i,...L2 , distinct within 

j = l...Mi,k=l...M2 

each of the four groups, with 

N + Li-Mi>N + L2-M2>0 (1.7) 

and such that the ?7's and /x's are not on the support curves of the bimeasure d^{x,y) 
in the x and y planes, respectively. (There is no loss of generality in assuming the first 
inequality in p.7p . since the results for the reversed case can be read off by symmetry in 
the two sets of variables.) The main result of this work is a new derivation of the following 
expression for the integral of a symmetric, rational function in two sets of variables 
{xa,ya}a=i,...N, having simple zeros at the points {^a}a=i...Li, {C/3}/3=i...L2 simple poles 
{Vj}j=i-:Mi, {/^fc}fc=i...M2 in terms of the biorthogonal polynomials {P„, S'„}, evaluated 
at the points {^a}a=i...Li, {C/3}/3=i...L2 5 bilinear combinations of these, together with 
the Hilbert transforms {Pn, Sn} evaluated at the points {fik}k=i...M2, {Vj}j=i-Mi- 

iSJ^ := I dfi{xi,yi) . . . I dfi{xN,yN)^N{x)AN{y) 

T^ niLi(e»-x.)nfci(C/3-^a) ^^^^ 

]V+L2-M2-1 JV+Li-Mx-1 



X 



e(Li, L2, M2, M2) W^/K WVh 



n=N n=N 



A^,(0A^,(C)AM,(r])AM2(/i) ' ' 



;i.9) 



where 



e(Li,L2,M2,M2) := (_l)i(Wi+M2)(Mi+M2-l)^_^^LiM2 ^^^^Q) 

and G is the (Li + M2) x (Li + M2) matrix 



G 



' N+L2-M2 N+L2-M2 

Kll {C,a,Vj) -^12 (^a,C/3) -P/V+L2-M2 (^a) ••• -P/V+Li-Mi-1 ('Co) 

N+L2-M2 N+L2~M2 

K21 ilJ'k,Vj) -^22 {^-kXp) -P/V+L2-M2(/^fc) ••• -P/V+Li-Afj-l (/^fe) 



:i.iii 



3 



with 

J ^'^ 

K,^{^,C):=Y,Pn{OSn{0 (1.12) 

n=0 

J -^-^ 1 
KlAi.V) ■■= Y^PniOSniv) + 7 (1-13) 

J -^-l 1 

K22{fi,C) ■■= J]Pn(/i)5n(C) + (1-14) 
n=0 ^ 

J -^-^ 

K^iifi^v) ■■= J2Pnif-i)Sniv) - Hifi,v) (1.15) 

n=0 

Similar expressions for the cases when one or both of the inequahties (jl.7|) is reversed will 
also be derived. (See Section 3.) 

Within the setting of two-matrix models, the integral (jl.9p may be interpreted as the 
expectation value of the product of Li evaluations of the characteristic polynomial of the 
first matrix and L2 of that of the second matrix, divided by similar products of Mi and 
M2 further evaluations of the respective characteristic polynomials. Integrals of this type 
were computed in the context of complex matrix models, where the variables {xa,ya) are 
replaced by pairs {za, Za) of complex conjugate values, in [T] for the special case when 
Ml = M2 = 0, and, more generally, in j2], for all values of {Li, L2, Mi, M2). The method 
of derivation used by these authors was based essentially upon recursive arguments, and 
is rather lengthy compared with the "direct method" that we present here. Other cases 
of the integral ()1.9|) . in which the condition ()1.7p does not hold, give rise to analogous 
determinantal expressions. These too will be derived in Section 3, using the same methods 
as those leading to ()1.9|) . 

Remark 1.1 Note that the multiplicative factor in front of the det^ term in ()1.9() is just 
a product of the inverse determinants that enter into rational interpolation formulae; e.g., if 
Li > Ml, L2 > M2, then 



(1.17) 



Al, (^) (C) Ami (??)Am2 (m) detGidetG2 
where 

Gi:= i £ i , ) l<a<Li, l<j<A/i , G2 ■■= ( Ca-p.,, Cfl) 1</3<L2, 1<J<M2 , (1-18) 

^ ' b<Li-l\Ii-l X P'fc ^' / c<L2-A/2-l 
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Remark 1.2 Expressions for the biorthogonal polynomials {p^^'''''^'^) (x), Sn'''''^'^^ (y)}nGN with 
respect to the modified measure 

together with their Hilbert transforms {Pn^''^'^'^\x), Sn''^'^'^\y)}ne'N may immediately be de- 
duced from formula ()1.9() . simply by making the respective replacements 

Li Li L2 L2 

n (^c. - Xa)^{x - Xa) JJ iia " Xa), JJ (C/3 " ya)^{y " ^a) ~ 

«=1 0=1 13=1 (3=1 

Ml Ml Af2 M2 

j=l j=l k=l k=l 



(2) 

in In and multiplying the result by the normalizing factor 



t(2) 
'■n + 1 



Remark 1.3 The corresponding formulae for the cases where one or more of the parameters 
within the groups {Ca}, {C/3}, {^j}' {t^k} coincide may be very easily determined from (|1.9jl . 
simply by taking the appropriate limits. This replaces the terms which appear multiply in the 
entries in (|1.1H) and the sums (|1.12|) - (|1.15|) by their derivatives with respect to the repeated 
parameters. 

As a "warm-up exercise", we also rederive the simpler analogous results for integrals 
of symmetric rational functions in one set of variables arising, e.g., in Hermitian one- 
matrix models. In this case, let {Pn{x)}n=o,i... denote the orthogonal polynomials with 
respect to a measure dfi{x), which may again, in general, be complex and supported on 
an arbitrary finite union of curve segments: 

dfi{x)Pn{x)Pm{x) = 6nm, (1-21) 



with leading term normalization 



P„(x) = ^ + 0(x"-^). (1.22) 



To guarantee their existence, the finiteness of the Hankel matrix of moments 

Mjk := J dfx{x)x^^^ < oo, 0, Vj, A; G N. (1.23) 
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is assumed, as well as the nonsingularity of all diagonal NxN submatrices {Mjk)o<j,k<N-i- 
The Hilbert transforms 

P4V) ■■= I d^^{x)^ (1.24) 
J 1]- X 

of the orthogonal polynomials are again assumed to exist. The partition function is 

N-l 

:= / d{fi)xi) ... dfi{xN)A%{x) = iV! n (1-25) 

The expression analogous to ()1.9|) for the integral of a symmetric, rational function in the 
variables {xa}a=i,...N, having simple zeros at the points {C,a}a=i...L, and simple poles at 
{rij}j=i...M, valid for N > M, is 

:= -L / df^ix,) ... [ df^ixM) n ^^#^f^A^(^) (1-26) 

^ a=l LLj=l\Vj ^a) 



" Ai(0AM(r7)nl~;i-MV^ V^iv-M(r7,-) ... P^+L-i(r7,)y'' 

where the two groups of complex parameters {^a}a=i,...L, {Vj}j=i,-M are assumed to have 
distinct values, the latter being evaluated off the contours of integration. The analogous 
formula for the same integral, valid for < M is 

J _ (-i)i^(^-^)(-i)^^^ v^nS-^ v^ ^.^ 

xdetf^f"! ••• p? . P \ (1-28) 

\PoiVj) ... PN+L-liVj) PoiVj) ■■■ PM-NiVj) J 

Such relations for polynomial integrands originated in the work of Heine and Christoffel 
fOlE] and were extended to the general rational case in ^H]. The case (jl.27p was recently 
rederived in the context of Hermitian matrix models by other methods [210 El . The direct 
method given here leads to relations (|1.27|) . (jl.28p in just a few lines. 

The key tool that is used in our "direct" approach is the following identity, which is 
just a multivariable partial fraction expansion for rational symmetric functions, valid if 
N> M: 



A^.(x)Am(^) ^ MN y y ^n-m{x[^]) 

lla=l lli=iv/j -^a) ^£5^^ ai<...<aM llj=ll'/o-j J-a^ J 
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Here x[a] denotes the sequence (xi, . . .xn) with the elements {xa^, ■ ■ -Xaj^j) omitted, and 

sgn((j) is the sign of the permutation a = ( ^ ' ' ' ^ ] G Sm- (To verify this identity, 

\ai ... (Tm / 

one simply notes that, viewed as rational functions in the ?7j's, the residues at all poles 
coincide, and both sides tend to as rjj^oo for any j.) Reversing the roles of {xa} and 
{rij}, an equivalent identity of slightly different form is valid, in this case for N < M: 

2^ sgn(a) ^ , (1.30) 



nil nSi(^. - ^a) (M - N)i ^ ^ uUv.. - ^a) 

Another relation that is of importance is the Cauchy-Binet identity from multilinear 
algebra. In invariant form, this states that if V is an oriented Euclidean vector space with 
volume form fi, and we have two sets of L vectors (P^, . . . , P"^), (S*^, . . . , S"^), then the 
scalar product of their exterior products (A^^^^P", A^^^S*^), defined by 

A^=i P" * A^S^ = (aLiP", a;^^,S^)Q (1.31) 

(where * denotes Hodge dual on exterior forms) is equal to the determinant of the matrix 
formed from the scalar products 

(aLiP", A^=i5^) = detG (1.32) 

In component form, if dim\^ = N + L, and the vectors {P", S°'}i<a<L have components 
{P", Sj^}i<j^k<L+N relative to a positively oriented orthonormal frame, this reads: 

,.l.......+......+.,n....fc^.+l.../c^.+.^l^^^ _ _ _ pj^^^^Sl^^ . . . St,, = iV!det(G), (1.34) 

where e denotes the Levi-Civita symbol and paired indices are again summed over. 

It is worthwhile noting that both ()1.29|) and p.31|) are particular forms of determinan- 
tal identities that may be deduced from Wick's theorem for products of free Fermi field 
operators. In a sequel to this work ^21, another method of deriving the relations ()1.9|) . 
fll.27|l and p.28j) is given, based directly on evaluation of vacuum state matrix elements 
of operators constructed from products and exponentials of Fermionic free fields. 
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2 The one matrix case. Proof of eqs. (11.271) . ( 11.281) 

To introduce the "direct" method used here, we begin by deriving the relation (jl.27j) . 
First, we recall the proof for the case M = given in 



^ / dfi{xi)... I rf/i(xAr) TT TT(^a - Xa)A^(a;) 

J J a=la=l 

dlJ,{xi)...J (i/i(xAr)A7v+L(x,0Aiv(x) 



1 Ln=N V " ' 
1 Ln=N V 



- j d^{xi)... j dn{xN)det{Pj{xa) Pj{^a)) det{Pkixh)) 



X P,,^,(ei) . . . P,,^,i^L)e'^ -'-PkAxi) ■ ■ ■ PkA^N). (2.1) 

Here the summation convention is used and eJ'i ---?'vJ'v+i ---?'v+i' and e^^'''^^ denote the Levi- 
Civita symbol in + L and N variables, respectively, the ranges of summation being 
< j, ji, • • • ,3n+l < N + L — 1 and < /c, fci, . . . , /c^ < — 1. Using the orthogonality 
relations ()1.2H) to evaluate the integrals yields 

liV = ^^-T — TTT det(PAr+a-l(^/3))l<a,/3<L- (2.2) 

To extend this result to the case of arbitrary M, we make use of the identity p.29|) . 
Substituting this into the integrand of (jl.27j) . using symmetry under permutations of the 
Xq's, and invar iance under relabeling of the integration variables gives 



ZNAMiv)iN - M)\ J r/i-zi J r]^ 



a=l j=l 



/ dn{Xi)... / dn{xN-M)^%-Mi^l^ ■ ■ ■ ^^N-m) Yl Yii^o, - ^a)Yl{Zj - Xa), 

a=l a=l j=l 

(2.3) 



where the combinatorial factor ^.^yj^j^^^, has been introduced to replace the sum over all 
permutations of ordered choices of M of the elements {xa}a=i,...,N, which are here relabeled 
{zi, . . . , zm), while the remaining ones are relabelled (xi, . . . , xn^m)- Now, applying the 
relation p.2|l with the L parameters ^ := (^i, . . . extended to L + M parameters 
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z) := (^1, . . . ,^L, zi, . . . , zm) to ()2.3p . and replaced by — M, and using ()1.25|) for 
both Zat and Ztv-m gives 



AL{O^M{v)Un=N-Mhn J Vl ' ^1 J Vj ' 

X det /^^^-^(^°) • • • ^N+L-ii^a) \ ^2 4) 

-M{Zj) ■ ■ ■ PN+L-l{Zj) J 

The multihnearity of the determinant when evaluating the integrals then gives (jl.27p . 

Now consider the case N < M. Substituting ()1.30p into the definition ()1.26|) of Ijv 
and using ()1.25p gives 

(_l)|A'(^-l) 

N\{M - N)\AL{O^M{r]) 

/ / d/i(xiv) — ^ 

y ^) 1 Ln=Q V 1 \.n=N V , 

m{M - N)\ Un=o VKA^iO^Miv) 

^/^v■■^/^AZi^--- 
det P,(r/.,^,)|o,.|M ^. . det ^^''^^ , (2.5) 



X 



-det ( Pj{rif^ ) o<j<M-N-i 



l<k<M-N 
<1<N+L-1 



X 



which is just the block matrix determinantal expansion of ()1.27|) . 



3 The coupled matrix case 

3.1 Case 1. N + Li - Mi > N + L2 - M2 > 

We now turn to the derivation of relation p.9p . which will be done along similar lines to 
the above. We again begin with the case with no denominator factors; i.e. Mi = M2 = 0, 
and Li > L^. For this case, using (jl.5p . 



'N a=l V"^ a=l (5=1 J 



9 



N-1 , \ A 

) A^+L,(x,OA7v+L2(z/5C) 



N 

X 

a=l ^ 

iV!A^,(e)A^,(C) 

X e^-^-'-^-'-^^^s^^iy,) . . . S,AyN)S,^,,{Ci) . . . ^fc.,,,(aj, (3.1) 

where in the first determinant < j, ji, . . . , Jat+Li < N + Li — 1, and in the second 
< k, k^+i, • • • , ^Ar+L2 — ^ + -^2 — 1- To complete the computation, we now make use of 
the Cauchy-Binet identity ()1.34|1 . To apply this to the expression ()3.1|1 . let L := Li and 
identify the vectors 5'^}i<a,/3<L as follows: 



pa 

j 



Pj-liU, l<^<Li, 0<j<iV + Li-l, 

Vi(C/3), 1<P<L2, 0<j <N + L2-1, 

0, 1<P<L2, N + L2<j <N + Li-1, 

Sn+m-ij, L2 + l<f3<L,, 0<j<iV + Li-l. (3.2) 



Using the Cauchy-Binet identity p.34|) and eq. p.Sp . the equality 1)3.11) gives the following 
expression for 

_ n=N V "-n 1 ln=Af V »-n ^ , /iV+L2 
"-N - (e)Ai (C) V (^a, 0) ^7V+L2(^a) • • • PN+L^-li^a 

(3.3) 

where 

N+L, ^±^-1 
^12(^,0:= 5^ PniOSniC). (3.4) 
n=0 

We now extend this result to the case of arbitrary (Li, L2, Mi, M2) satisfying ()1.7p . 
We detail the derivation only in the case when the stronger inequality 

N > max(Mi, M2) (3.5) 

holds. For the intermediate cases, when N lies between Mi — Li and Mi or between 
M2 — L2 and M2, formula ()1.9|) may be derived by similar computations. 



10 



Substituting the identity p.29|) for both denominator factors Y[a=i 11^=1 ('7i ~ ^a) and 
rifciil/^fc - Va) into the integral (HIHI) defining I^^ gives: 



I(^) = -A^— --- J]sgn(a) 5:sgn(a) (-1)^-^"^ E (-l)""^^^'^ 




(3.6) 



In this sum we must distinguish: 

m := the number of a^-'s that coincide with b^s 

Ml — m = the number of aj's that do not coincide with any fe^'s 

M2 — m = the number of b^s that do not coincide with any a/s 
Note that, if Mi + M2<N,m can vary from to min(Mi, M2), but if < Mi + M2, it 
can only take values m > Mi + M2 — N. The number of distinct ways in which two such 
ordered sets ai < . . . < aMi, and 61 < . . . < 6mi with exactly m common elements can be 
chosen from the numbers (1, . . . , A^) is: 



In view of the invariance of the integrand in (jl.9p under permutations of the pairs 
{{xa,ya)}a=i,...N, and the freedom to relabel the integration variables, we may express 
the integral as 



iN 

m,Mi,M2 ■~ 



{N - M1-M2 + m)!(Mi - m)!(M2 - m)!m! ' 



(3.7) 



^ \ i (Ml +M2 ) (Ml + Afa - 1) 



min(Mi,M2) 



l(2) 



Z]^^AMi(r7)AM2(yu) 





m=max(0,Mi+M2-Af) <7&Sm-^ 



X 



X C, 



X 




AM2-m(^) 



W 
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>^ n 

a=l 



a=l k=l 



M2-m L2 Mi-m 

Yi (^fc ~ n*^^^ ~ n ~y<^) ] ^ 



(3=1 



^ ^iV-i\/i-Af2+m(^)^iV-i\/i-i\/2+m(l/) (3-^ 

where we have made the following changes of notation in the integration variables 



1/1, 

•^m+l ) ■ 
Vm+l, ■ 



• 5 Xjri 
■ ? l/m 

; 3; Ml 



?/Mi+l, • • • , UKh+hh-m 
XMi+M2-m+lj • • • ) 3; AT 
2/Afi+M2-m+l ■ ■ ■ ,yN 



Zl, 
Wi, 

Zl, 
Wl, 

fl, 

Wl, 
Xi, . . 

yi ■ ■ 



) Zm 
, Wm 

) ZM2—m 

, WM2-m 
, XN-Mi-M2+m 
, yN-Mi-M2+m 



(3.9) 



In determining the sign factor in the first line of ()3.8p , we have replaced the sums Y2j=i % 
and J2k=i appearing in ()3.6|) by their values for the case (oi, . . . , a^/J = (1, . . . , Mi) and 
(61, ... , 6^2) = (1; • • • ) ""^^ Ml + 1) • • • ; Ml + M2). This leaves a residual factor (—1)™^ = 
(— 1)™", giving an alternating sign in the sum over m. 

A further simplification can be made in ()3.8|) by noting that, in the sums over the 
elements of the symmetric groups 5^/^ and Sm27 all terms in the integrations over the 
{zi, Wi}i=i,„m, {zj,Wj}j=i,„Mi-m and {zj, Wk}k=i...M2-m Variables coming from pairs of per- 
mutations (a, cr) for which the sets {{ai, ai)}i=i^,„m, {crm+j}k=i,...Ah-m and {^Mi+k}k=i,...M2-m 
are invariant contribute the same value to the sum, and there are m!(Mi — m)\{M2 — m)\ 
of these. That is, all left cosets of the subgroup Sm,MiM2 '■= SMi-m x SM2-m C Sm^ x 
5*4/2 that permute separately the first m elements in both (1, . . . , Mi) and (1, . . . , M2), 
the last Ml — m elements in (1, . . . , Mi) and the last M2 — m elements in (1, . . . , M2) con- 
tribute the same term in the sum. Hence we may choose one representative [cr, a] G {Smx x 
SM2) / Sm,Mi,M2 from each coset, multiplying this term by the factor m!(Mi— m)!(M2— m)!. 

We may also now apply relation ()3.3j) . with the replacements N^N — Mi — M2 + m, 
Li-^Li + M2 — m, ^2-^/^2 + Ml — m, as well as the expression ()1.5|) for both Z^'' and 



r(2) 

^N-M-i-M2+m 



to evaluate the integrals over the {(x^, ya)}a= 



■■l...N-Mi-M2+m 



variables. Using 



the resulting sum becomes 

(_l)|(J^'A+Af2){Afi+A/2-l) Y^N+L2-M2-1 ^^^Yl 



l(2) 



N+Li-Mi-1 
n=N 



iV!Ai,(0Az.2(C)AMi(r/)A 



min(Afi,Af2) 

X E ( 

m=0 



E 



sgn((Tjsgn((T 




dfi{zi,Wi 



Utli^o 



Ut lie, 



Wi 
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X det 



X det 



N + L2-M2 _ _ N+L2-M2 ~ ~ ~ 

Ki2 {Zk,Wj) K12 (ZkXfi) PN+L2-M2{^k) ■■■ PN+Li-hh-l{Zk) 

N+L2-M2 _ N+L2-M2 

K12 {C,a,Wj) K12 (^a,C/3) PN+L2-M2{ia) ■■■ PN+Li~Mi~l{^a) 

( ^\hiMi+M2){Mi+M2-l) T-fN+L2-M2-l ^ /JT J-\^+^^~^'^^-^ . /TT 
y lln=Af V ""n 1 Lm=N V "-n 

Ai,(OAi,(C)AM,(r/)AM,(/i) 

min(Mi,M2) m 

m=0 [a-,cr] i=l 

■N+L2-M2 N+L2-M2 

^21 if^aMj+k^^'^m+j) ^22 (/^5Mi+fc'C/3) 'PjV+L2-M2(/^^Mi+J • • • 'PiV+Li-Mi-l(/iaM^+J 
N+L2-M2 N+L2-M2 

^11 i^a,V(7^+j) K12 {iaXp) PN+Li-A-hi^a) ■■■ PN+Li-Mi-li^a) 



(3.10) 



where 



^n(e,r]) := / df,{x,y)&^^^^^ K,2{^,y) (3.11) 
/C22(/^, C) := / df^ix, y) ''^-''^' K,2ix, C) (3.12) 
K2\{[i,r]) := / d^{x,w) / dfi{z,y) Ki2{z,w) 

(3.13) 

n{^^ V) := j d^x, y) (^_^)(^_^) (3-14) 

Vnifi) ■■= [ d^{x,y) ^'^~^^^^ — —Pn{x). (3.15) 



Here, in the first hne of eq. (j3.1(J|) we have cancelled the combinatorial factor m!(Mi — 



m)!(M2 — m)\ with the corresponding expression in the denominator of C^^^ »^ and 



used the relation (jl.Sp for 'Z^^N-Mi-h-h+m ^'^ cancel the further factor (A^ — Mi — M2 + n)!, 
leaving only the normalization factor H^if Y{n=NX^L^-~M2 numerator. 

Note also that, although the range of summation in m is from to min(Mi, M2), in the 
case when N < Mi + M2, all terms with < m < Mi + M2 — N vanish, because the 
determinant factors, which are of dimension (Li + M2 — m) x [Li + M2 — m), have entries 
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X det 



that are formed from scalar products of vectors of dimension N + Li — Mi < Li + M2 — m, 
and hence have less than maximal rank. 

We now note that the sum over cosets [a, o"] G (5*^1 x SM2) / S(m,Mi,M2) fl^.lUp is an 
expansion, as a sum over homogeneous polynomials in the terms {7i(/ijt, 77,)} i<j<Mi , of 

l<k<M2 

the single (Li + M2) x (Li + M2) determinant 

('_lU(Mi+M2)(Mi+M2-l) T-rA^+L2-M2-l fT i-rAr+Li-Mi-l rr- 
j(2) _ I ^) lln=jV y ""n. iin=N V 

Ai,(0Ai,(C)AM,(r/)AM2(/i) 
' N+L2-M2 N+L2-M2 

^21 il^kiVj) ~ T^il^k^Vj) ^22 il^kXls) VN+L2-M2{l^k) ■ ■ ■ VN+Li-A4i{l^k ) 
N+L2-M2 N+L2-M2 

A^ll {ia-,Vj) K12 iiaXp) PN+L2-M2{ia) ■■■ PN+Li-Mii^a) 
( -\\UMi+M2)(Mi+M2-1)( A\LiM2Y['^+L2-M2-l 1 i-rAf+Li-Mi-l rr— 
_ \ I '-J lln=Af '''n 1 Ln=N+L2-M2 V 

Az.,(OAz.2(C)AM,(r7)AM2(/^) 

(Ar+L2-M2 N+L2-M2 

A^ll {C.a,Vj) -^12 (^a,C/3) PN+L2-M2{^a) ■■■ PN+Li~Mi{C.a) 

N+L2-M2 N+L2~M2 

^21 {fik,r]j) -'H{fik,r]j) ^22 ifJ'kXls) VN+Li-Miif^k) ■ ■ ■ VN+Li-Ah{f^k) 

(3.16) 

Notice also that, by adding linear combinations of the last Li — Mi — L2 + M2 columns 

N+L2-M2 N+L2-M2 
of the matrix in (j3.16|) the terms { /Cn {C,a,i]j)^ ^21 {f^k,Vj)} be replaced by 

Af+Li-Afi iV+Li-Afi 

{ ^11 (^a,^?i), /C21 (/ifc,?7j)} to give 

/ -|\i(A//i+Af2)(Afi+Af2-l) rr^+^2-Af2-l /T- T-rAT+Li-A/i-l rr- 
t(2) ^ I -"-J^ lln=Af V ''■ra 1 ln=N V 

Ai,(OAi,(C)AM,(r7)AM2(/^) 

(N+Li-Mi N+L2-M2 

A^ll {C.a,Vj) K12 {iaXp) PN+L2-M2{ia) ■■■ Pn+Li-MiH 

N+Li-Mi N+L2-M2 

^21 {Hk,Vj) -'Hil^k,Vj) ^22 ifJ'kXls) VN+L2~M2il^k) ■ ■ ■ ^Af+Li-Afi (/ifcj 

(3.17) 

As a final step we note that by separating the integrands in (j3.1H) - (j3.15j) into the sum 
of their principal parts at the poles x = t] and y = fi and polynomial parts of degrees 

Af+L2-Af2 

< Li — 1 in X and < L2 — 1 in and using biorthogonality (which implies that K12 
is the kernel of an integral operator projecting onto the first N + L2 — M2 biorthogonal 
polynomials), the integrals (prTTD - (nrTH|) . for J = + Li - Mi or J = + L2 - M2 may 
be reduced, at the values {^a, C/?? Vj^ f^k}, to the following 

N+Li-Mi 111 N+L^-M-i 



^11 = YYiia- Vj) Kii {^a,Vj) (3-18) 



0=1 
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N+L2-M2 il N+L2-M2 

/C22 (/^fe,C/3) = ||(C/3 - Atfc) K22 (/ifc,C/3) (3.19) 
/3=1 

AT+Li-Mi I± ^ N+Li-Mi 

(/ifc, rij) - nink, Vj) = Vj) ii(C/3 - yUfc) ^21 (/ifc, r]j) (3.20) 

0=1 /3=1 

Vnifik) = ll{(p- f^k)Pn{fik), n > L2. (3.21) 

Substituting these expressions into the determinant term in ()3.17|) . factoring out the 
diagonal matrices diag(l, . . . , 1, Y\%i{Ci5 - /^i), • • • , Y[%i{Cp - l^Ah)), making the replace- 

N+Li-Mi N+Li-Mi N+L2~M2 N+L2-M2 

ments Ku {^a,r]j), ^21 ilJ^k,r]j)^ (^a,'7i), ^21 (/Ufc, 77^) in the first L2 + 

Ml columns of the determinant, by adding linear combinations of the last Li — Mi — L2 + 
M2 columns and factoring our the diagonal matrices diag(]^^^^(,^Q, — r/i), . . . naLi('^a ~ 
r^Af J, 1 . . . , 1)) on the left and right when evaluating the determinant then gives the rela- 
tion (fOU) . 

Remark 3.1 Although the computation was done here for the case when N > Mi,N > M2, 
the other cases, in which Mi — Li < N < Mi, or M2 — L2 < N < M2, or both, may be derived 
similarly, leading to the same formula H1.9() . 

Formulae analogous to (IT^ for the cases + Li - Mi > 0, + L2 - M2 < and 
+ Li — Ml < 0, + L2 — M2 < may similarly be deduced using the second form of 
the partial fraction identity ()1.30p whenever the inequalities relating the degrees of the 
numerator and denominator polynomials in ()3.6|) require it. These are derived in the next 
two subsections. 

3.2 Case 2. N + Li - Mi > > N + L2 - M2 

In this case the integral I]y is given by the following determinantal expression: 

y-rN-Mi+Li-l /T — y-r M2- N-L2-I IT — 

= 6(Li, L2, Ml, M2)^^ -^^^ 

" A.,(OA..(C)AM.(r^)AM.(/.) ^'-''^ 

where 

e{Li L2 Ml M2) "= — l)f^(^-l)+l^i(^i+l)+l^2(i2-l)+MiAf+M2Li+MiL2+LiL2 ^3 23) 
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and G is the (Li + M2) x (Li + M2) matrix 



Pb{^o) 



G ■= det ( '^^^^^ _ j ('3 24) 

\Hif^k,Vj) ^m(/ifc) y ' 

Here the row indices are, sequentially, 1 < « < Li and 1 < k < M2, and the column indices 
1 < j < Ml, 1 < /3 < L2, < 6 < iV - Ml + Li - 1 and < m < M2 - iV - L2 - 1. Note 
that the first three column blocks of ()3.24|) coincide with those of the matrix G defined 
in p.llll if one understands the orthogonal polynomials Pn{x), Sn{y) and their Hilbert 
transforms to vanish for negative n. Formula ()3.22|) is valid whenever N + Li — Mi > 0, 
N + L2 — M2 < 0, but it is easier to demonstrate assuming the stronger conditions 

A^>Mi, A^ + L2-M2<0, (3.25) 

which is what we do in the following. The intermediate case, when Mi — Li < N < Mi 
may be proved through a similar computation. 
Expressing I^"* in this case as: 

du(x V ) UtUi^a - ^a) 

X (3.26) 

and applying identity ()1.29|1 with respect to the (x, 77) variables and ()1.30|) with respect 
to the ((?/,C)),Ai) variables gives 

^ iV!(M2 -N- L2)!(nl'o' K)AlAO^mAv)^mM 




X 



J2 (-1)^^=1 ^ sgn(a) ^ sgn(a 



'^M2-L2-N{VaN+L2 + l^ ■ ■ ■ '''^o-Afa) 



ai<-<aMi fTg5A/i ?e5M2 11/3=1 (/^trjv+fi C/?) 

n 



^ dfl{Xa,ya) utilise, - Xa) \ I:^N -M^A^^]) 



\a=\ 



,_^YI,N+\{N+L2){N+L2-l)+\Kh(M,+l) ]-[L2^ HSll/^fc " 0) 



(M2 - iV - L2)!Mi!(iV - Mi)!(nl-; K)/\i^,{0/^l2{0^mAv)^aM 



X 
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X 

a=l 



&]^^L^^k^ 0<m<M2-N-L2-l, l<k<M2-N-L2 



b 

Q, / l<a<N-Mi, l<a<Li 
0<b<N-Mi+Li-l 



L2 (3-27) 

^_^Y4,N+\{N+L2){N+L2-l)+\M,(M,+l) 1 (/"fc " Cp) 

{M2-N- L2)m\{N - Mi)!(nl~o' hn)ALA0^L2{0^MAv)^AM 

^■'—^ J J/ I 1 l<a<JV-A/i, l<ci<Li 

""SS'a/j^ u&Sm2 0<i)<JV-Afi+Li-l 

n ( TT^ _ /■ I ^^^^^^^N+L2+k>0<m<M2-N-L2-l. l<k<M2-N-L2^ (3.28) 

o_i \Po-JV+;3 S/3/ 



X 

/3=1 



where 



J {r]-x){fi-y) 
f 

Xbifi) := J dfi{x,y) ^_ ^ . (3.30) 



We now note that ()3.28p is just the block determinant expansion of 
I 



(nl"o /^n)Az.,(OAi,(C)A,,,(r^)AM,(^) 



ln=0 

xdet ^'''f'* "" l (3.31) 



0<6<iV-A/]^-Ll-l, 0<m<M2-iV-L2-l 
l<a<Li, 1<I3<L2, l<j<Mi, l<fc<A/2 



(The combinatorial factors (Ms - L - L2)!Mi!(iV - Mi)! in (jT^ are cancelled by the 
multiplicity with which the sub determinant factors occur in the sums over a G Smi, 
a G Sm2-) a further simplification can be achieved by applying elementary column 
operations. To see this, we take the integrations appearing in the matrix elements of 
fj3.3ip outside the determinant, re- writing it as: 



I 



(nlV ^n)Az.,(0AL2(C)AM,(r])AM2(/i) 



ln=0 

xnf/^^V^tf%fe^ 4 ^ .-y (3.32) 



fc=l ^ \ " "Sa 
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Now, viewing as a rational function in Xk with a simple pole at Xk = Tjj, we 

may re-express it as the sum of the pole term plus a polynomial of degree < Li — 1 



riaLll^a ~ ^k) _ YlgLli^a ~ Vj) , \^ a 7 /o qq\ 



where 

Li-l Li 



E ^^7C= 11(^/5 -^^■)- (3.34) 



7=0 /3=1 



Since N — Mi +Li — 1 > Li — 1, all monomials in Xk of degree < Li — 1 appear in the second 
column block of the determinant in ()3.32|) . We may therefore add linear combinations of 
these columns to those in the first block to obtain the equivalent expression 



det I fe-O M.-C^ = det Z"^-"*) f 

e / I -n/^i(e/3-r],) e 



nn&-''.)d'=M V! f' V ) 

JV-Mi+Li-l M2-JV-L2-I ii M2 / 1 p /™ \ 1 c /,, ^ 

nv^ n v^nn&-''.)det(v! v 't' 

n=0 n=0 a=l j=l ^ Vj-^c t>\^aj 



(3.35) 

where further elementary column operations were made in the last lign to replace the 
monomials x^,^^ and /i^ by the corresponding biorthogonal polynomials. Substituting 
this expression into ()3.32|) and evaluating the integrals gives 



I 



(.) _ (-l)^^^^+l(^+^^)(^+^^-^)+I^MM.+l) ^^^^ ^M.^(^^ _ ^^.) ^ ^M^^^^^ _ 

X U/K. nV^det^ ^^C. ^ J. (3.36) 

n=0 71=0 Vj Ka ' 

Finally, reordering the rows and columns suitably, bringing the third column block into 
the second position and interchanging the two row blocks, we arrive at the expression 

dm. 
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3.3 Case 3. iV + Li - Mi < 0, iV + L2 - M2 < 

(2) 

In this case the integral Ijy is given by the following determinantal expression: 

^ n=Q V't'nl [n=0 V '^n 

^ mil nail - u^, u^Mc, - /^.) , . . 

where G is the (Mi + M2 - iV) x (Mi + M2 - iV) matrix 



G := det 



PiiVj) I , (3.38) 

V ^ ^ I 




with rows labelled consecutively by 1 < a < Li, 1 < k < M2, and < £ < Mi — Li — N 
and the columns by 1 < j < Mi, 1 < /3 < L2 and < m < M2 - L2 - iV. 
To derive this formula, we begin by expressing the integral in the form 



ly = — ^ 7^^ / dfi{xi,yi) . . . / d^{xN,yN) 



,(2) ^ u:uiru^a-v,)UT=iU, Mc, - /^.) 

z£)a^,(OAz.,(C) 

^ Ajv+L,(x,OAjv+L,(l/,r/) 

nil n'Li n;i(e« - - ^a) u^, uttic, - -vaY 

Using the identity (jl.3(J|) twice, first with the Xa variables replaced by the combined set 
(xq, ^j) and N^N + Li, and also for the set {i/a, (k) with N—*N + L2, we obtain 



j(2) 



^ iV!(Mi - Li - iV)!(M2 -L2- N)\{Yf^:^ K) Ai,(0Ai,(C)AM,(r7)A„,(/i) 

^ ^ ^ AAfi-Li-Af(?7CTAr+i^+l, • • • ,'7fTAfi)AAf2-L2-Af(/^5jV+L2+15 • • • ^/^SMa) 



TV 

n 



dfii^Xay Ua 



(_l)L,(iV+MO+L2(7V+M2)+|L.(L.+l)+iL2(L2+l) ^^^^ ^M.^(^^ _ ^^^^ jjM,^^^^ _ 

N\{M, - Li - iV)!(M2 - L2 - iV)!(nl'o ^n) A^,(0A^2(C)AM,(r/)A„2(/i) 

aG5A/i 5G5A/-2 \a=l / \a=l'l''f^+" \/3=l ^'^'v+zs <./3^ 
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(_l)L.(A.+MO+L.(A.+A/.)+iL,(L,+l)+iL.(L.+l) J^-L^^ ^^^^^ _ ^L^^ ^M.^^^^ _ 



(nl"o ^n) AL,(OAL,(C)AM,(r/)A^.(/x) 



N 




xdet ^ 1, (3.40) 
V r/J 

where the rows are labelled consecutively by 1 < A; < M2, 1 < a < Li and < i < 
Ml- Li-N and the columns by 1 < j < Mi, 1 < /3 < L2 and < m < M2 - L2 - A^. 
Finally, interchanging the two row blocks, and replacing the monomial entries rf^ and /i^ 
by the biorthogonal polynomials Pj{rjj)k) and Smil^k) respectively, with suitably modified 

————— ('2') 

normalization factors, we arrive at the expression (HOTIl for V^'. 
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